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Self-SynchronizingAdaptive Multistage Receiver for Wireless 
Communication Systems 

CROSS-REFERENCE TO RELATED APPLICATIONS 
5 This patent application claims the benefit of the filing date of United States 

Provisional Patent Application SeriarNo. 60/248,469, filed November 14, 2000 and entitled 
"Adaptive Multistage Filter-Bank Receiver for Wireless Communication Systems," the entire 
contents of which are hereby expressly incorporated by reference. 

10 FIELD OF THE INVENTION 

The present invention relates to a filtering system and method for wireless 
communication signals. More specifically, the present invention is directed to a method and 
system for adaptively filtering asynchronous signals in wireless commimication systems. 



1 5 BACKGROUND OF THE INVENTION 

Direct-sequence code-division multiple-access (DS-CDMA) has become an important 
component of recent communication systems. This is due to its many attractive properties for 
wireless communications such as the one depicted in FIG. lA. In DS-CDMA 
communications, users' signals are multiplexed by distinct code waveforms, rather than by 

20 orthogonal frequency bands in frequency-division multiple access (FDMA), or by orthogonal 
time slots in time-division multiple access (TDMA). All users can transmit at the same time 
with each of them being allocated the entire available, broad-band, frequency spectrum for 
transmission. This method of channel allocation alleviates the problems commonly 
associated with other forms of channel allocation such as unused resources. 

25 In CDMA, demodulation requires the suppression of two forms of noise, ambient 

channel noise, which is modeled as additive white Gaussian noise, plus multiple access 
interference (MAI), which is highly structured. The conventional receiver, which uses a 
simple correlation of the received signal and its user-spreading sequence (a standard matched 
filter), does not take into account the existence of MAI. The MAI caused by any one user is 

30 generally small, but as the number of interferers and their total power increases, the effects of 
MAI can become substantial. These effects fall within, what is called, the near-far problem. 

The conventional receiver is severely limited by the effects of such MAI. One way to 
alleviate the MAI problem is to use power control schemes. However, power control is not 
easy to accomplish, especially in a wireless environment, where the power levels often vary 
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dramatically. Even if perfect power control is used, the system may still need to be 
implemented with a forward error-correcting code to obtain a good performance. Another 
way of handling the near-far problem is to use near-far resistant receivers, which, in principle, 
can yield a substantial improvement over the standard receiver with perfect power control. 
5 Multiuser detection techniques were proposed many years ago to .mitigate the MAI 

problem. For example, S. Verdu, in "Minimum probability of error for asynchronous 
gaussian multiple-access channels," IEEE Trans, on Inform, Theory, vol. IT-32, pp. 85-96, 
January 1986, [1], the entire contents of which are hereby expressly incorporated by reference, 
proposed and found the optimal multiuser detection criterion and maximum-likelihood 
10 sequence detection for DS-CDMA. Unfortunately, Verdu 's detection and demodulation 
methods were generally too complex to realize in a practical DS-CDMA system. 

As a consequence, a number of simplified detector-demodulators have been proposed 
to obtain near optimal performance with reduced computational complexity. (For example, 
see R. Lupas and S. Verdu, "Linear multiuser detectors for synchronous code-division 
15 multiple-access channels," IEEE Trans, On Inform, Theory , vol. 35 , pp. 123-136, 
January 1989 [2]; R. Lupas and S. Verdu, "Near-far resistance of multiuser detectors in 
asynchronous channels," IEEE Trans, on Commun., vol. 38, pp.496-508, April 1990 [3]; Z. 
Xie, R. T. Short, and C. K. Rushforth, "A family of suboptimum detectors for coherent 
multiuser communications," IEEE J. Select. Areas Commun, , vol. 8 , pp. 683-690, May 
20 1990 [4]; U. Madhow and M. L. Honig, "MMSE interference suppression for direct-sequence 
spread spectrum CDMA," IEEE Trans. On Commun., vol. 42, pp. 3178-31888, De- 
cember 1994 [5]; M. L. Honig and J. S. Goldstein, "Adaptive reduced-rank residual 
correlation algorithms for DS-CDMA interference suppression," in Conference Record of 
the Thirty-Second Asilomar Conference on Signals , Systems & Computers , 1998 , vol. 1, 
25 pp. 551-555 [6]; and D. J. Lee, Adaptive Detection of DS/CDMA Signals Using Reduced- 
Rank Multistage Wiener Filter, Ph.D. thesis, University of Southern California, March 2000 
[7], the contents of which are hereby expressly incorporated by reference.) 

These receivers, however, deal primarily only with detectors that require a precise 
knowledge of the propagation delays of all users, which is usually unknown a priori. To use 
30 such an algorithm, the propagation delays need to be estimated and as a consequence, the 
receivers suffer from errors that occur with the estimation of the propagation delays. (See S. 
Parkvall, E.G. Strom, and B. Ottersten, 'The impact of timing errors on the performance of 
linear ds-cdma receivers," J, Select, Areas Commun,, pp. 1660-1668, October 1996. 
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[8], the contents of which are hereby expressly incorporated by reference). To obtain a better 
performance, the transmitters and receivers needed to be synchronized. Moreover, the 
transmitters and receivers, in this case, need to be synchronized in order to obtain a better 
performance. However, the synchronization requirement adds more complexity and cost to 
5 the system. 

The conventional acquisition system that utilizes a tapped delay line with its weights 
being the elements of the spreading sequence also suffers from the same near- far problems. 
Many near-far resistant timing-acquisition algorithms have been developed. One approach is 
to jointly estimate the needed parameters of all users, as shown in Y. Steinberg and H. V. 

10 Poor, "On sequential delay estimation in wide-band digital communication systems," in 
Proc, Intl. Symp. Information Theory, San Antonio, TX, 1993, p. 423 [9]; and S. Y. Miller 
and S. C. Schwartz, "Parameter estimation for asynchronous multiuser communications," in 
Proc. Conf. Information Sciences and Systems, Baltimore. MD, 1989, pp. 294-299 [10], 
the contents of which are hereby expressly incorporated by reference. While such techniques 

1 5 can provide excellent results, they are computationally complex. 

Recently, maximum-likelihood synchronization of a single user is developed in S. E. 
Bensley and B. Aazhang. "Maximum-Likelihood synchronization of a single user for Code- 
Division Multiple-Access communication systems," IEEE Trans, on Commun., pp. 392-399, 
March 1998 [11], the contents of which are hereby expressly incorporated by reference. 

20 However, the method in [11] requires a training period. This requirement is not feasible in 
certain situations. For example, whenever there is a significant change in the MAI level, 
during a deep fade or a tuming-on of some near-in strong interferers, retraining is required, 
and the data transmission need to be interrupted accordingly. Similar limitations apply also to 
the minimum mean-square-error (MMSE) interference suppression techniques for the timing- 

25 acquisition, developed in U. Madhow, "MMSE interference suppression for timing 
acquisition and demodulation in Direct-Sequence CDMA systems," IEEE Trans, on 
Commun,, pp. 1065-1075, August 1998 [12], the contents of which are hereby expressly 
incorporated by reference. 

By the use of a subspace decomposition technique, such as multiple signal 

30 classification (MUSIC), another type of timing estimator is developed in E. G.Strom, S. 
Parkvall, S. L. Miller, and B. E. Ottersten, "Propagation delay estimation in asynchronous 
direct-sequence code-division multiple access system," IEEE Trans, on Commun., Vol. 44, 
pp. 84-93, January 1996 [13]; and S. E. Bensley and B. Aazhang, "Subspace-based channel 
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estimation for code division multiple access communication systems," IEEE Trans, on 
Commun., vol. 44 , pp. 1009-1020, August 1996 [14], the contents of which are hereby 
expressly incorporated by reference. This type of estimator does not need a training sequence 
and provides good results when unknown information sequences are sent. However, this 
5 timing estimator still involves extensive computations due to the need for an eigen 
decomposition. Also, this technique assumes that the number of active users is fixed and 
needs to be known or estimated. As a result, there is a need for a system and method for a 
relatively simple, accurate, and cost effective filtering of an asynchronous wireless signal. 

10 SUMMARY OF THE INVENTION 

An adaptive near-far resistant receiver for wireless communication systems, such as, 
asynchronous DS-CDMA systems is developed without the use of prior synchronization. The 
only requirement is knowledge of the spreading code of the desired user. Also, there is no 
need of a training period or of data-free observations. 

15 In one aspect the present system describes a method and system for filtering of an 

asynchronous wireless signal comprising the steps of receiving a data vector; of using the 
received data vector to update the weight coefficients of an adaptive filter without prior 
knowledge of synchronization and without a need for a separate training period or data firee 
observations. 

20 In another aspect, the present system describes an adaptive near-far resistant receiver for 

an asynchronous wireless system comprising: means for receiving an asynchronous data 
vector including a spreading code; using the received asynchronous data vector, means for 
updating weight coefficients of an adaptive filter without prior knowledge of synchronization 
of the spreading code of the data vector; using the updated weight coefficients, means for 

25 determining synchronization of the spreading code; and means for demodulating the output of 
the filter using the determined synchronization of the spreading code for obtaining a filtered 
data vector. 

In yet another aspect, the present system discloses a method performed by a computer for 
filtering interference and noise of an asynchronous wireless signal comprising: receiving an 
30 asynchronous data vector including a spreading code; using the received asynchronous data 
vector, updating weight coefficients of an adaptive filter without prior knowledge of 
synchronization of the spreading code; using the updated weight coefficients information to 
determine synchronization of the spreading code; and demodulating the output of the filter 
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using the determined synchronization of the spreading code for obtaining a filtered data 
vector. 

The method further comprises the step of dividing the data vector represented by x[/] 
into two channels Xi[/] and dx[i] using a transformation Ti on x[/], represented by Tix[/], 
5 wherein the transformed data vector x[/] does not contain information about a designated 
sender's spreading code 5i, and d\[i] contains primarily only information about the spreading 
code Si and residual data from correlation of 5| and x[/]. 

In yet another aspect, the present system describes an adaptive receiver for filtering 
interference and noise of an asynchronous wireless signal comprising: means for receiving 
10 an asynchronous data vector including information data bits; means for updating weight 
coefficients of an adaptive filter without prior knowledge of synchronization of the 
information data bits for determining the start of the information data bits; and means for 
demodulating the output of the adaptive filter. 

1 5 BRIEF DESCRIPTION OF DRAWINGS 

Other embodiments of the present invention will become readily apparent to those 
skilled in the art from the following detailed description. As will be reahzed, the invention is 
capable of other and different embodiments and its several details are capable of modification 
in various respects, all without departing from the spirit and scope of the present invention. 
20 Accordingly, the drawings and detailed description are toi)e regarded as illustrative in nature 
and not as restrictive. 

FIG. 1 A is an exemplary block diagram of a mobile phone channel; 

FIG. 1 is an exemplary block diagram of an asynchronous DS-CDMA system; 

FIGs. 2(a)-2(g) are exemplary timing diagrams of construction of a typical DS signal; 
25 FIG. 3 is an exemplary block diagram of a receiver according to one embodiment of the 

present invention; 

FIG. 4 is an exemplary block diagram of a structure of a test statistic; 

FIG. 5 is an exemplary block diagram of a multistage decomposition of a likelihood ratio 
test for M=4; 

30 FIG. 6 is an exemplary decomposition of data vectors of FIG. 7; 

FIG. 7 is an exemplary conceptual block diagram of a self-synchronized detector; 
FIG. 8 is an exemplary flow diagram for deriving test statistics by the iterative version of 
the multistage adaptive filter; 

5 
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FIG. 9 is an exemplary conceptual flow diagram of block 77 of FIG. 7; 
FIG. 9 A is an exemplary algorithm for finding the maximum of |Re{y[/]}| over NS 
period; 

FIG. 10 is an exemplary flow diagram for demodulation system of FIG. 7; and 
5 FIGs. 1 1-20 are performance charts for various embodiments of the present invention. 

DETAIL DESCRIPTION 

In this invention, an adaptive near-far resistant receiver for asynchronous DS-CDMA is 
described which assumes no prior knowledge of synchronization. The receiver of the present 

10 invention only needs to have knowledge of the spreading code of the desired user. Also there 
is no need of a separate training period or of data-free observations. The receiver is 
composed of a near-far resistant acquisition criterion and demodulation scheme. 
Furthermore, a lower complexity version of the receiver is developed that is related to the 
concept of the multistage Wiener filter described in J. S. Goldstein, I. S. Reed, and L. L. 

15 Scharf, "A multistage representation of the Wiener filter based on orthogonal projections, " 
IEEE Trans. Inform. Theory, vol. 44, pp. 2943-2959, November 1998 [16], the contents of 
which are hereby expressly incorporated by reference. This results in a substantial reduction 
of the computational burden and a rapid convergence in an adaptive implementation. Finally, 
an iterative version of the multistage decomposition is developed. 

20 As shown in L. E. Brennan and I. S. Reed,"A digital adaptive anteima for JTIDS," The 

final report on contract F19628-88-C-0098 , Adaptive Sensors, Inc., December 1989, 
submitted to the Electronic Systems Division AFSC, USAF [15], the contents of which are 
hereby expressly incorporated by reference, acquisition is an optimum signal detection 
problem analogous to radar-signal detection. In [15], the detector of the signal looks at each 

25 sampling interval for a new data set of N samples that occur periodically, where N is the 
spreading gain of the system. This is accomplished by a moving window or its physical 
realization, using a tapped delay line (TDL). This moving window is used to sweep through 
each periodic time interval of N samples. In each such time period, an estimate of the timing 
is obtained from the maximum of a set of optimum adaptive filter outputs. Thus, 

30 synchronization can be done sequentially over the code length by the maximum likelihood 
detector. For demodulation, an adaptive MMSE receiver described in[4] and [7] is 
implemented. One approach for implementing the receiver is the equivalent of a TDL, finite 
impulse response (FIR) adaptive filter bank. Such an FIR detector-demodulator adaptive 
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filter bank is converted then to the considerably simpler bank of iterative or infinite impulse 
response (IIR) filters. 

The authors in [16] have developed a new type of reduced-rank adaptive Wiener filter. 
With the concept of the reduced-rank multistage- Wiener adaptive filter a lower complexity 
5 version of the receiver is developed. As a consequence, covariance-matrix inversions or 
eigen decompositions are avoided, thereby making possible a substantial reduction in 
computational complexity. The multistage scheme not only eliminates matrix inversions, but 
its adaptive version which is developed in [7] can be realized directly from the raw A/D 
converted digital data without any need for an explicit estimation of the covariance matrix. 

10 These important features contribute to a significant reduction of the computational cost and 
the amount of sample support usually required to accurately estimate the covariance matrix, 
as shown in L S. Reed, J. D. Mallet, and L. E. Brannan, "Rapid convergence rate in adaptive 
arrays," IEEE Trans, Aerosp. Electron. SysL, AES-10, pp.853-863, November 1974 [24], the 
contents of which are hereby expressly incorporated by reference. 

15 In this invention, a new iterative method or an IIR version of the reduced-rank 

multistage- Wiener adaptive filter is developed. In this new method, the coefficients of each 
filter of a filter bank are updated every time that a data vector is received. Thus, each filter 
continuously keeps track of the changes of the environment. Furthermore, this new algorithm 
is implemented more efficiently because it obviates the necessity of matrix multiplications. 

20 This new algorithm can be realized as an iterative or IIR adaptive filter bank with each filter 
of the type suggested in [6], whereas the original algorithm in [16] or the training-based 
algorithm in [7] needs to be realized as an FIR adaptive filter. 

For simplicity reasons, only binary phase-shift keying (BPSK) modulation is developed 
and used in this disclosure, however, it will be recognized by those skilled in the art that other 

25 modulation schemes for other embodiments of the invention are also possible. Therefore, the 
receiver can be applied to many other modulation schemes that utilize other symbol alphabets 
and spreading-sequence generators. Furthermore, the present invention is not limited by the 
type of modulation code or the code length of the desired user or of the other users. For 
example, the family of 4-phase sequence developed in P. V. Kumar, H. F. Francis Lu, T. 

30 Helleseth and D.-J. Shin, "On die Large Family of Low Correlation Quaternary Sequences 
S(2;. " Proceedings of the IEEE Intnl Conf Personal Wireless Comm., (IGPWCIOOO), 
Hyderbad, pp.33-37, December 17-20, 2000 [26], the contents of which are hereby expressly 
incorporated by reference, can also be used as a spreading sequence in addition to the Gold 
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code (described in Appendix A). Synchronization and demodulation of the receiver system is 
derived in detail in the following sections. The reduced-rank version of the receiver that is 
based on the multistage representation of the Wiener adaptive filter is also described below. 

5 Asynchronous DS-CDMA System Model 

Consider an asynchronous DS-CDMA system that uses BPSK modulation with K 
simultaneous users on an additive white Gaussian noise (AWGN) baseband channel, as 
illustrated in FIG. 1 . The received signal can be written in the form, 

10 r(0= ^^^kMs,{t-mT,-r,) + n{t), (1) 

m--oa k = \ 

where the parameters and other quantities are defined as follows: 

Tb Information symbol interval. 
15 n{t) Additive white Gaussian noise (AWGN). 

bk[m] The w-th information symbol of the ^-th user and bk[ni] € {-1,0,1} 

where symbol zero denotes the fact that the ^-th user is off (has either not started 
or has finished transmission). 
Ak Amplitude of the ^-th user. 
20 Tit . The propagation delay of the ^-th user. 

Here, Sk{t) is the spreading waveform of the ^-th user, given by 

Sk(t) - I; cf<l>{t-jT,), (2) 

where Tc is the chip interval and <j3{t) is the chip waveform, which is assumed to approximate 
25 a rectangular pulse of duration Tc. (Each bit in the spreading code is called a chip.) The 
spreading code waveforms have unit amplitude and are assumed to be periodic with the time 
period Tf,, In one symbol period, there are N = TJTc chips, modulated with the spreading 

code ( \ . . . , c^^.j ). Here, also N is called the spreading gain. 

To understand how the asynchronous DS-CDMA signal is determined, an example for K - 
30 2 and = 31, is presented next. In the DS-CDMA system, the information-bearing sequence 

8 
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for each user is spread over a broader bandwidth by means of a spreading wave-form that is 
unique to that user. FIGs 2(a)-2(g) depict construction of a typical DS signal. FIGs. 2(a) and 
2(b) are parts of two information sequences, FIGs. 2(c) and 2(d) are two spreading 
waveforms, FIGs. 2(e) and 2(f) are the two modulated waveforms to be transmitted, one for 
5 the desired user 1 and the other for the interfering user (user 2). FIG. 2(g) is the combined 

waveform of FIGs. 2(e) and 2(f). Note that user 2 is assumed to be delayed by t^-t^ ~ 4T^ 
seconds with respect to user 1 . FIG. 2(a) shows two bits of an information sequence for user 
1. The spreading waveform of user 1 is shown in FIG. 2(c). It is constructed by modulating 
the chip waveform ^/) by means of a spreading code. In this example, a Gold code of length 

10 31 is employed for the spreading code of user 1, given by Si = (-1, -1, -1, 1, 1,-1, 1, 1, 1, 1, 
-1, 1, 1,-1, 1,-1,-1,-1, 1, 1, 1, 1, 1, 1,-1, 1,-1,-1,-1,-1 ,-1) (see Appendix A). Each 
bit in the spreading code Si is called a chip. FIG. 2(e) shows the modulated baseband DS 
signal that corresponds to the multiplication (modulation) of the information sequence and the 
spreading waveform in FIGs. 2(a) and 2(c). The interfering information signal of user 2 that 

15 contributes to the received signal, is assumed to be offset in time from that of user 1 by 4Tc 
and is illustrated in FIG. 2(b). Also, FIGs. 2(d) and (f) show the spreading waveform and the 
modulated signal of user 2. The spreading code for user 2 is given by 52 = (1, 1, 1, 1, -1, 1, 
-1, -1, -1, 1, -1, -1, 1-1,-1, 1,-1, -1, 1, 1, -1, 1, -1, -1, -1, -1, -1, 1, -1, -1, -1). The 
baseband transmitted signal for the time interval [t, /+2rj, shown in FIG. 2(g), is the sum of 

20 the modulated signals of both users. 

To facilitate the derivation of the receiver, the communication channel is assumed to 
approximate a time-invariant channel over some interval of time exceeding N chips, i.e. the 
channel is quasi-stationary. In particular, the number, the powers, the delay profiles of the 
users and the interferers are assumed to be fixed (but unknown). Only the turning on-and-off 

25 . of the information sequences, their b^s, and the AWGN are changing in time. In scenarios of 
this type, it turns out that the receiver is able to adapt quickly to any such changing 
environment. As a consequence such assumptions affect the usefulness of the receiver of the 
present invention only minimally. 

An exemplary receiver structure is shown in FIG. 3. For simplicity in its realization, the 

30 proposed receiver is described using a baseband-equivalent structure. Such a baseband 
complex signal process is physically achieved by the combination of quadrature demodulation 
and a phase-locked loop (PLL). (For example, see J. G. Proakis, Digital Communications, 
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McGraw-Hill, New York, NY, 3"* edition, 1995, Chapter 6 [25], the contents of which are 
hereby expressly incorporated by reference). This converts the received radio frequency (RF) 
modulated signal to a baseband complex-valued signal. This means that after converting the 
received signal to baseband complex- valued signals, the received signal is passed through a 
5 filter that is matched to the square-wave chip (or other type of waveform that is used to be the 
chip waveform). The output of this chip matched filter is 

x(t) = l^^' r(tW = f r(t'-t)dt'. (3) 

10 After that the output of the chip matched filter is sampled once every 7; seconds, where S = 
TJTs is an integer and S > 1 . This means, that the sampling rate is greater than or equal to the 
chip rate. Then the output is used as the input to the AT-element tapped delay line (TDL). The 
signal is sampled at time zT, so that the tapped delay line can be expressed by the column 
vector, 

15 

X W = [x(iTs). x{iTs - Tc ) , ... , x(iTs - (N - l)Tc)] , (4) 

where [.]^ denotes matrix transpose. Similarly the tap-weights' vector is expressed as the 
column vector, 

20 

w \i] = [ [i], [t\, . . W] ^ • (5) 

The components of the weight vector w[/] are determined later in Eq. (30) of Section 3.1. 
The output of the tapped delay line filter, y[i], is given by the inner product of the vectors 
25 in Eqs.(4) and (5) as follows : 

y\i] = w [iV x[i-]= X w„[i\'x(iTs-Tc + niTt)r (6) 

where [.]* denotes the conjugate (Hermitian) transpose of a matrix, and [.]* denotes the 

30 conjugate of a complex number. This output is passed through the time acquisition and 

10 
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demodulation system to obtain the information bit after synchronization. This time 
acquisition system can be modeled as a filter bank constructed of NS parallel receivers, each, 
of the type shown in FIG. 3, in order to identify the time phases at which to demodulate the 
information symbols of the desired user. 

Receiver Structure. Synchronization System 

As shown in [15], the synchronization problem in communications is a signal 
detection problem analogous to radar detection. This problem can be modeled as a binary 
hypothesis problem, where hypothesis Hi corresponds to target signal presence, and 
hypothesis Hq corresponds to target-signal absence. Thus at each time sample of the A/^ vector 
x[i] the synchronization receiver need to distinguish between two hypotheses of some user, 
say user 1. The target signal under hypothesis Hi is denoted by the A^-vector ^i^iSi, where 

Ai is the amplitude of user 1, bi is the information bit of user 1, and Si = [Cq\c[^^, . . . , 

^N-i y is signature sequence of user 1, where the observation vector at the sample time /T, 
is the N-vector x[/], given in Eq. (4). 

The noise environment is composed of A/- vectors that are due to the combined effects 
of the interference sources and white noise. The noise process n(t) is assumed to approximate 
complex Gaussian noise with zero mean by an application of the central limit theorem [11] 
and a covariance matrix, defined as follows: 



Thus, at a particular symbol interval the two hypotheses, which the adaptive detector 
need to distinguish at each sampling time, are given by 



given b\ where b\ = +1 or -1. Note that the condition bi = 0, x[i] is equivalent to the null 
hypothesis //q. 

For simplicity, let N = Tt/Tc and S = TJTs , where both numbers are integers. If the 
spreading code modulation is approximately band-limited, the case S = 2 would correspond 
with Nyquist sampling. It is easy to show that the covariance matrix of the sample vector x[/] 




Ho 



X [i\ = n[z] , 

X [i] = AibiSi + n[i] . 
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is periodic with the integer time period NS, In other words, x[i] is a cyclostationary process 
with period NS, 

The random vector x[/], when conditioned on the first user's information symbol bu is 
approximately Gaussian under both hypotheses. The condition probability densities of x[/] 
5 are 

(8) 

,-x'li]R„M'M>"l 
;r^det(R„[i])' 



10 

where det(R) denotes the determinant of R. From Eq. (8), the conditional probability density 
of p(x [i] I Hijbi) can be expressed in terms of p(x [i] | Hi,bi) for 6, = ±1 as follows: 

P(x [«] I Hi) = Xp(^i 1 ^1 >i) ' (^) 

15 

where p(6, = 1) = p(6i = -1) = 1/2. Then, the likelihood-ratio test (LRT) takes the form, 

A = PWillilil > y, (10) 
P(x[«]|Ho) < 

where y is the detection threshold. The substitution of the conditional probability density 
functions, given in Eq.(9), into Eq.(lO) yields: 

20 



2 



pix[i]\H,,b, =l) + p(x[i]|H,,b, =-1) 
P(x[«l|Ho) 



g-(i^,5i)'R„ '['l(>^i»i) [g-2Re((/\,5,)'R„ ^ g2Rel(/iiSi)'R„-^[i]x[i]) | 



(11) 



Thus, the likelihood ratio test in Eq.(lO) reduces to 

12 
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Hi 

A = cosh (2Re{(A,b,s,) * R„ » [«]x[i]}) J Yi , (12) 

Ho 

where yi is a new decision constant. Since the hyperbolic cosine function cosh(.) is a 
monotonically increasing function in the magnitude of its argument, the test in Eq.(12) is 
clearly equivalent to 

H, 

\Re{q[i]}\ = |Re^|R„-'[,]x[,]] > 72, (13a) 



where 



9[0 = s|R„-'[«]x[«], (13b) 



and 72 is another constant. 

The test statistic in Eq. (13a) can be used at each clock time or phase within the time 
10 period NTc to test for the existence of the signal of a particular user, If 

Q[i\ = [ |Re{t/[i-NS + I]}|} (14) 

is the set of such tests in increasing time sequence of the form, given in Eq. (13a) over a past 

15 spreading code period total of NS clock times, then the integer clock time i e {/, z-1, . . . , 
A^5+l}, such that 

Re{fl[i]}| - max I Re{(7[i - A:]} I (15) 

is attained, is the time that synchronization is most likely to occur within the interval (z- 
A^5+l, /). The "greatest of 'sequence-in-time test in Eq. (15) detects sequentially the 

20 maximum of all likelihood tests in the filter bank of the tests in the set Q [i] in Eq. (14). 

The matrix Rn[0 in Eqs. (13a) and (13b) is the covariance matrix of the signal-free 
observation AT- vector x[/] = n[/]. In the classical adaptive-array application it is normally 
possible to collect signal-free observation vectors. However, in a CDMA-type system, 
obtaining signal-free observation vectors at the receiver is not at all a trivial task. The test in 

25 Eq. (13a) needs to be considerably modified in order to be realized with both signal-plus- 
noise and noise-only data vectors. 



13 
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To accomplish this modification test in Eq. (13a) is expressed equivalently next in, what 
is called, the generalized sidelobe canceler (GSC) form. Consider an N x N nonsingular 
matrix Ti with the structure. 



Ti = 



(16) 



where B, is, what is called, a blocking matrix whose rows are composed of any orthonormal 



basis set of the nullspace of the normahzed signal vector Ui = s,/-y/slsi . That is, 
10 Biui = Bisi = 0. (17) 

The test statistic in Eq. (13a) can be re-expressed in terms of the transformation T, , defined in 
Eqs. (16) and (17) as the expression. 



15 



Re{slR„-' [i]x[i]} = |Re{s;T;(T,R„[ilT;)-'T,x[«]} 



(18) 



To evaluate the latter test statistic in Eq. (13a), let the transformed output vector be expressed 
in the form, 



20 



Tix[il = 



'n\x{i] 




'd,[i] 


B,x[i] 







(19) 



where [/] = u*x[i] is a signal-plus-noise scalar process and x,[/] = B^x[^], is an (AT- 1) - 
dimensional process with the signal removed. 

Here, the data vector x[z] is divided by the transformation Ti into two channels or 
25 paths, namely, and xi[q, respectively, as shown in FIG 4. The transformed data vector 
Xi[/].does not contain any information about the signal Si because it was obtained by a 
projection of the data onto the orthogonal complement of the signal subspace. Actually, in 
the channel of the dx[i] process, obtained when the data vector x[i] is projected onto the signal 
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10 



subspace, the interference has been suppressed already due to the low cross-correlation 
property of the other users' spreading codes. This dx[i] process is the same process that is 
obtained from the cross-correlation conventional detector which yields for low interference 
level an acceptable performance by the exploitation of this property. However, with severe 
MAI due to the near-far problem, it is possible that the contribution of the interference in the 
di[i] channel is still significant. Thus, the other channel Xi[z] which contains information only 
about the interference-plus-noise can be used to partially cancel the interference-plus-noise 
process in the di[i] channel. Then the herein described multistage decomposition for a 
Wiener filter provides a natural and optimal way to accomplish such a step-by-step 
cancellation process. 

The correlation matrix of the transformed vector Tix[/] can be expressed as the following 
partitioned matrix: 



15 



20 



25 



T,R,[/]T; = 



u;RJz]u, u1rJz]b1 
_B,RJz]u, B,RJz]Bl 



<[0 CM] 



where, 



alii] = ulRJiK = Edrfjif}, 
,^,,,[1] = B,Rji]u, = E{xJiK[0}, 



(22) 



(20) 



(21) 



(23) 



with E{.} denoting the expected value operator. Next the signal-free correlation matrix R„[/], 
needed in Eq.( 1 3a), can be epressed in terms of R,[j] by 



RnW = RxW - A?s,s; . 



(23a) 



where s,s| is the x A^matri outer product of vector Sj by itself. If one defines the positive 
number, = ^sj Sj , it can be shown, using Eq.(23a), that 
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T,R„[i1T[ = 



(23b) 



With the aid of the matrix inversion lemma in Jerry M. Mendel, Lessons in 
Estimation Theory for Signal Processing, Communications, and Control, Prentice-Hall, 
Englewood Cliffs, New Jersey, 1995 [17], the contents of which are hereby expressly 

incorporated by reference, the inverse of TiR„[/]T| is computed to be 



10 



(t,r„[ot;)-^ 



= (^i[i] + A?<5f)-' 

(24) 
where 



1 <,Ji]r;V[i] 



15 



Then, the transformed signal TiSi is given by the vector. 



(25) 



Tis, = [V(s[sJ 0...0]* =S^[1 0 ...0]* , 



(27) 



20 Hence, the test statistic in Eq. ( 1 3a) can be re-expressed as 



Re{s|R;'[f]x[i]}| = |Re{s;T;(T,RJilT»)-^T,x[f]}| 

= |Rek [i\ + AlS^rs;in\ - wj^ [i]B, )x[f ]j , (28) 



where. 



25 



(29) 
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From Eq.(28) the term 5| is a positive constant which can be moved to the other side of the 
test in Eq.(13a). Also the term i^^i]-^ A^S^)'^ defined in Eq.(26), does not vary 
significantly over the period NS when the signal AySi is comparable to the interference-plus- 
noise power level so that [0 + ^.^5,^)'* is approximately a positive constant over the 

5 same NS time period. Hence also the term [i] + (J,^ can be moved to the other side 
of the test in Eq.(13a). Thus the test statistic q[i] in Eq.(13a) can be modified to a new test 
statistic given as follows: 

Nymil = |Re{(ul-w^[z]B,)x[z]}| J ya, (30) 

Ho 

where the weight vector in Eq. (5) has been replaced by 
10 u;-w^[/]B,, (30a) 

and y[i ] which replaces q[z] in Eqs.(14) and (15) is defined by the modified likelihood ratio 
test 

y[i] = (u;-w^[z]B,)x[z], (30b) 
The block-diagram structure of this new test statistic is shovm in FIG. 4. Evidently, this 
15 test statistic has the form of the classical sidelobe canceler, used in the past to suppress or 
cancel interferers or jamming of radars and communication systems. Also note that the 

modified test to determine i for synchronization in Eq.(15) is given by 



20 



|Re{y[i1}| = max |Re{y[z-fc]}l (30c) 

with q[i] replaced by y[i] in Eq.(30b). 

From the test statistic in Eq.(30) the expression wgsc is the only term to utilize the 
covariance matrix. It is easily verifiable that 



25 R,,[i] = B,R„[/]B| = B,(Rji]-A?s,s;)Bl = B,RJi]B; ,(31) 

= B,R„[/]s, = B,(Rji]-A'saS;)Si = B,R,[i]s, . (32) 



These relations in Eqs. (31) and (32) demonstrate that the test statistic in Eq. (30) remains 
invariant to the presence or absence of a signal component in the observation vector x[il. 
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Demodulation System 

Once synchronization is established, data demodulation is accomplished by the use of 

a MMSE detector [4, 5]. Let x[i] , where i is the time phase or clock time defined in Eq. 
(15) at which synchronization occurs, be the observation vector obtained upon 
5 synchronization. Then, the estimate of the information symbol bi is obtained fi-om the 
estimation statistic in Eq. (30) has the form, 

b, = sgn(Re{(A,s,)*R;^[f]x[?l}). (33) 

where sgn denotes the sign operator. Since the ampHtude Ai and the norm of the signal 

sequence, namely -^s* is positive, the estimate of 6i can be simplified to 

10 b, = sgn(Re{u;R;^[Ox[f]}). (34) 

However, for the modified test statistic in FIG. 4, the estimate of bi is obtained by replacing 
the weight vector u* R;^ [i ] in Eq.(34) by the new weight vector in Eq.(30). Thus 

15 b, = sgn(Re{(u;-w^[i]B;)x[f]}). (35) 

is flae modified estimate of the information symbol to use hereafter in this invention. It should 
be noted that the term (u* -wJ^[i]B Jx[i] is used in common with both the 
synchronization system and the demodulation system. This term can be computed and stored 
20 during the acquisition process. It does not need to be recomputed for demodulation. 

An exemplary implementation of the Receiver .-^ 

Recall that x^'"'[i] is the observation vector at the sampling time iTs of the m-th 
symbol. Assume that L approximately independent snapshots (samples) of the observation 
25 vector x^'"*[i] are available for the initial acquisition of the detector parameters. These data 
are given in matrix form by 

Xo[«]A [x<>)m,...,x(^'[f]]. 
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The maximum likelihood estimator for ['] is computed then to be 

Kin = jix^-W-ni] = yXol'KW- 
The two statistics needed to find wcsccan be computed by the following: 

5 r,,,, =B,RJ/]s, =iB,Xo[/]X*[/]s, 

By substitution of these two estimates into Eqs.(30) and (35), an adaptive version of the 
proposed receiver is obtained. 

Lower Complexity with the Multistage Representation 

10 Recently, rank-reduction methods have received considerable attention in the context 

of a Wiener filter [16] and in J. S. Goldstein and L S. Reed, "Reduced-rank adaptive 
filtering," IEEE Trans. Signal Processing, vol. 45, pp. 492-496, February 1997 [18], the 
contents of which are hereby expressly incorporated by reference, in which form the detection 
tests in Eqs.(30) and (35) can be applied. To achieve this, the observation vector is projected 

15 onto a lower dimensional subspace, and the Wiener filter is constrained to lie in this subspace. 
This procedure makes possible optimal signal detection and accurate signal estimation while 
allowing for a lower computational complexity and a smaller sample support. 

In this section a reduced-rank version of the receiver is developed which is based on a 
new Wiener filter structure, called "multistage representation of the Wiener filter," [16]. This 

20 multistage structure represents a pioneering breakthrough in that it achieves near full rank 
detection and estimation performance with often only a small number of stages. This can 
result in a substantial reduction of the computational burden and a rapid convergence in an 
adaptive implementation without any need for matrix inversion. 

To achieve the desired multistage decomposition of the likelihood ratio test in Eq. (13a), 

25 successive orthogonal decompositions are applied to the observed data. Actually, the 
derivation of the test in Eq, (30) fi*om Eq. (13a) can be considered to be the first stage of the 
multistage decomposition. To accomplish the second stage of the decomposition, the same 
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technique is applied to the term Xi[/] = B|X[i] and wgsc defined in Eq.(29). First, define U2[i] 
as the vector. 



U2[l] A 



(36) 



10 



where, r^^^^ [i] is the correlation defined in Eq. (22). Then, define a new nonsingular trans- 
formation T2[i], similar to Ti in (16), by 



T2W = 



(37) 



where the rows of the (N-l) x matrix B2H is the null space of U2[i], i.e., B2[i]u2[i] = 0. 
Thus, the new transformed vector under transformation T2[i] is represented by 



15 



20 









BJilxJi], 







(38) 



in accordance with the definition in Eq. (19) for the first stage decorhposition. The procedure 
for the decomposition of the likelihood ratio test in the first stage is repeated in the second 
stage. By following the same arguments that lead to Eqs. (28) and (29), the second stage 
decomposition of the likelihood ratio is calculated fi-om Eq. (28) by the relations 

y[il = (u;-w[3sc[ilBi)x[i] 

= (ul - oy; [i](u» [t] - mK I'l^a I'DBi M'l ) 
=dJi]-ty;[i](dJil-4j,]R;^Ji]xJ,]), (39) 



25 where (o,[i] A JJ^I , = ^^.^ [i]^,,, [i] and ^2[i1 = < \i] - r^^,^ [iIRJ [«K,., W • Also, 
define an error signal €3 ['] ^ 
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e^li] = d,[i]-ri^,^mlli]x,[i]. (40) 
The variance of £2 [i] in Eq. (40) is given by 

5 

= E{e,\i]el[i]}. (41) 

Thus, the substitution of Eq. (40) into Eq. (41) yields: 

10 <[/] = E{(djf]-4ji]R;^ji]xjz])(dji]-4joR;U'>2[n)*} 

= ai [i\ -2 mK \^K^. + C.i'i^x, [']Rx, UK ['K,., [«i 

= < [^1 - ['IRx, [iK,, \i] = m . (42) 

Recall from Eq. (25) that = al ['K' UKd^ I'l = < t'l - ^•^^'"^ I'l = 

15 rfj [j] - i^^a^ [i]R'l [i]x, [i] is defined similarly to [i] in Eq. (40). From Eqs. (39) and (40), 
one obtains the relation 

ejf] = y[«] = eji] . (43) 

20 between e, [i] and €2 [i] . Hence the variance in Eq. (25) can be expressed further by 

m = c^', [']T^(T2R,,[t]T»)-^T2i;.,,, [I] 

= C7l\i]-^;'[i]S^\i] , (44) 

25 thereby directly relating the variance with the corresponding variance ^2[i] of the second 
stage of the multistage decomposition. 

A continuation of this decomposition process, extending Eq. (39), yields the TV-stage 
likelihood ratio test in terms of only scalar quantities in recursive form as follows: 
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y[t\ = d, \i] - co\ [z \i\ - co\ ... \i]{du-x \i\ ' [i]d^ 1^]) - )) • 



(45) 



More complete details of the multistage decomposition are developed in [16]; and J. S. 
5 Goldstein, 1. S. Reed, and P. A. Zulch, "Multistage partially adaptive STAP CFAR detection 
algorithm , " IEEE Trans. Aerosp.^ Electron. Syst., vol. 35, pp. 645- 661, April 1999 [19], 
the contents of which are hereby expressly incorporated by reference. 

Rank reduction is realized as a truncation of the full-scale structure of the formula in Eq. 
(45) at the M-th stage, where M<N. That is, the rank-M detector is obtained by stopping the 

10 decomposition at stage M , i.e. by setting Bj^i] 0. As consequence , [i] = and 



The complete recursion procedure for the rank-M version of the likelihood ratio test in 
Eq. (13a) is summarized in Table 1 as a pseudocode (also given in [16] and FIG. 5). The 
exemplary block diagram of FIG. 5 illustrates the likelihood ratio test for M=4, 



^mW = <['•] = < m- 



15 



Table 1 shows recursion formulas for the reduced-rank likelihood ratio test. 
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Table 1 

Initialization 



10 



15 



20 



ui = , B, = nuil(u,) , Xo[j] = x[q 



Forward Recursion 

For j = 1 to M - 1 



dj[i] = u;.[/l x^..Ji] 



X;[i] = B.[«]x..Jz] 



null(u^.^Ji]) 



End 



Define x^_,[i] 
Backward Recursion 

l|2 



25 



<W = E{||d^[«]|| } 
For j - M- 1 to 2 

End 



= E{||rf,.[il|| ) 
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Adaptive Realization of The Multistage Decomposition of The Test 

Usually, the covariance matrix for an adaptive filter is not known a priori with any 
5 certainty. Thus to perform the detection test in Eq. (30), it is necessary to find a statistical 

estimate R^[z] to substitute for R^[i] . In addition, the inversion of R^[i] is needed for each 
zTs where 0 < z < NS - 1 . This involves high computational complexity, and it is not a 
trivial task to accomplish. 

However, the multistage structure that is introduced in the previous section not only 

10 eliminates the task of covariance matrix inversion, but its adaptive versions can be realized 
directly firom the observation data without the explicit estimate of the covariance matrix, 
given in Eq. (47). Furthermore, when the code vector of the desired user is known to the 
demodulator, the coefficients of the adaptive multistage detector, i.e., Uj, By and fij, are 
determined without any knowledge of the reference information symbols. 

15 In [6, 7], a training-based algorithm for the multistage decomposition is introduced. 

Using this algorithm, one does not need to calculate an estimate of the covariance matrix. 
The coefficients of the detector can be estimated directly fi-om the received data vectors and 
then be used to demodulate the data after the training period. Also, in this algorithm, the 

dimension of the blocking matrix can be kept the same (i.e., in an AT x matrix form) 

20 for every stage. Evidently, such a blocking matrix can be obtained fi-om the following 
formula: 

B^n = i-u,.[i]u;[«], (48) 

where Uy[z] is on estimate of imit vector Uy[/] in Table. 1. 
25 In this manner, the lengths of the registers needed to store the blocking matrices and 

vectors can be kept the same for every stage, a fact that is very desirable for either a hardware 
or software realization. To obtain this algorithm let 

u;mX^._J«], (49) 
B,[tlX._Ji]. (50) 
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d][i] A [df^\i],...,df^[i]] 
X,[i] A [xf[i],...,xf'[i]] 



24 



42121/RRT/R535 

The estimate of the cross-correlation vector i^ .d i^] computed then as 




Also let the estimate variance of dj [i] be computed by 

^ m=l 

Table2 depicts a training-based multistage decomposition algorithm. 
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Table 2 : 

Assume that L independent samples of the observation vector x^'"^[/] are 
for the initial acquisition. Let [i] A [ x^^^ [i ] , . . . , x^^ [z ] ] 
Forward Recursion 

Initialization : Let vl^ - and = I - xx-u^- 

For 7=1 to(Af-l) 



Calculate d- and x^. 



X,.[i] A [xf[i],...,xf>[f]] = B.[z]X^_Jf], 



Calculate the (j + l)-th stage basis vector. 




Determine the (/ + l)-th blocking matrix B 



End 



Compute rf)^' [i] and set it equal to [i] 



Backward Recursion 



For 7 = Af - 1 to 2 



Estimate the variance of 
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End 



Estimate the variance of 



Calculate the y-th scalar Wiener filter 
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Then the variance of the error can be obtained from the difference equation, 

5 The rank-M algorithm shown in [7] is summarized in Table. 2. For the blocking matrix 
B^.[i], given in Eq. (48) in Table. 2, one observes that the matrix multiplications required to 
compute X^ [i] are given in Eq. (50). These matrix multiplications can be avoided by the 
substitution of into Eq. (50). That is, 

X.[z] = B^mX._Jz] = (I-i,[i]i;W)X^._Jz] 

10 = - a X -_Jz] = X..Jf] - 

Using this result, a further simplified version of Table 2 is given in Table 3. Note that this 
new structure no longer requires the calculation of a blocking matrix and the computational 
complexity is reduced significantly. 

In this invention, an iterative algorithm of the multistage decomposition is developed. 

15 Instead of estimating the parameters in the batch mode, the parameters also can be estimated 
iteratively. By this means, the coefficients of the filter are updated every time that data are 
received. Thus, the filter can more rapidly keep track adaptively of the changes in the 
interference environment. 

20 An Iterative Version of the Multistage Adaptive Filter 

In order to afford the possibility of following the statistical variations of the observed data 
when the detector operates in a non-stationary environment, a weighting factor is often used 
to ensure that data in the distant past are forgotten. A special form of weighting that is 
commonly used is the exponential weighting factor. As shown in Simon Haykin, Adaptive 

25 Filter Theory, Prentice-Hall, Upper Saddle River, New Jersey, 1996, the contents of which 
are hereby expressly incorporated by reference [20], an estimate of the cross-correlation 

vector i; .J . [z] can be expressed as 

=Z(i-«)"-"'xf>[f]d;.'">[i]*. (53) 
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where a is a positive constant less than 1 . The inverse of a is, roughly speaking, a measure 
of the memory of the algorithm. The special case, when a approaches zero, corresponds to 
infinite memory. 

Table 3 shows a simplified training-based batch-mode multistage decomposition 
algorithm. 
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Table 3: 

Assume that L independent samples of the observation vector x^'"-' [i] are available 
for the initial acquisition. Let X^[i\ A [ x^^^[i] , . . x^^*[z] ] 
Forward Recursion 

Initialization : Let = and x^li] = x[/] 

For j = 1 to (M - 1) 

Calculate d- and x^ 

dj\i\ = u;[t]x^.J«] 

djM A [df\i],..., df>[r]] = X^_,\i\, 

A [x;^)[i], x'/'m] = x,._,[i] - i^m djw 

Calculate the (/ + l)-th stage basis vector, 
1 



u..Ji] 
End 

Compute d^M^[i\ and set it equal to e';;" [0 
Backward Recursion 

^ m=l 

For j = M- 1 to 2 

Estimate the variance of d-{i\ 



30 



42121/RRT/R535 



*J m=1 



'm=l 

Estimate the variance of e 



J 

2 r,i £-1 



Calculate the y-th scalar Wiener filter by 0)j[i] 



End 
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The summation on the right hand side of Eq. (53) can be written as the difference or iterative 
equation. 



= il-ayri';p]+xf[i]df^[ih (54) 
5 where r^"^"^^ [i] is the previous value of r^"] . [z] , and the term x j"^ [i]df^ [i] is obtained from 
the new data. Hence, for each iteration, only two quantities (r^"J^^[z] and x*."*[z]dj"^[i] ) 
are necessary for updating r^"] . [i] . The length of the observation-data that contributes to the 

estimate, r^"] [z] , is determined by the weighting factor a A ^ , where L is the time 

constant corresponding to the length of the batch-mode ip Tables 2 and 3, 
10 The exponential weighting factor is substituted in the algorithm in Table 3. This results 

in the iterative algorithm shown in Table 4. The algorithms in Tables 3 and 4 do not require 

any need for matrix multiplications. Also in order to find the estimate i , the time phase 
interval in Eq. (30c) at which the information of user 1 is found, the algorithm in Table 4 is 
reahzed as an IIR adaptive filter bank while the algorithms in Tables. 1, 2 and 3 utilize FIR 
15 adaptive filter banks which need substantially more computer memory and computation to 
realize. 

The time phase acquisition process to find the estimate i , defined in Eq. (30c), is 
illustrated by the conceptual block diagram shown in the FIG. 7. The composition of the data 
vectors in the exemplary block diagram of FIG. 7 is shown in FIG. 6. The observation data 
20 vectors are obtained by sliding a window of length through the original data. The test 

statistic x/^''^[i]can be derived by the subroutine shown in the FIG. 8. In addition, the 

conceptual block diagram to find the value of / = i that achieves the maximum value of 
|Re{y[/]}| is found in FIG. 9. An exemplary algorithm for finding the maximization in 
Eq.(30c) is shown in FIG. 9A. Finally, FIG. 10 shows the block diagram that is used to 
25 demodulate the signal of the desired user. 
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Numerical Results 

In this section simulations were conducted to demonstrate the performance of the pro- 
posed asynchronous CDMA receiver. An asynchronous BPSK DS-CDMA system with 
the number of users = 6 is considered in these simulations. The spreading sequence of each 
user is a Gold sequence of length A^ = 3 1 (see Appendix A). A single desired user, say user 
l,is acquired and demodulated in the presence of multiple-access interference. Table 4 
depicts an exemplary pseudo-code of the iterative multistage decortiposition algorithm. 
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Table 4: 

Let k be the it-th clock time, where A: = 1 is the first time the data is observed. Also 
let a be the memory or time constant of the algorithm. 
For ^ = 1 to 

5 = si, ui'^[i] = = hi, x^o'^[i] is the received vector. 

Forward Recursion 

Fory = 1 to A/-1 

<'[«•] = u;.*^W x<!> [z] 

1 0 Calculate the (/' + 1 ) - th stage basis vector, 

rji 

End 

Backward Recursion 

For 7 =Mto 2 

Estimate the variance of [i] 
20 Calculate the j - scalar Wiener filter Gff\i\ by 

Calculate the (7 - l)-th error signal e^*] [i] by 

End 
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The output at time fc-th is y^*' [i] = e'/* [i] 

End 
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The power ratio between each of the five interfering users and the desired user is 
designated the MAI level. In this simulation all interfering users are assumed to have the 
same power ratio relative to the desired user. This power ratio is denoted by a quantity called 
the near- far ratio (NFR), defined by 



where the subscript 'k' denotes user A^, A: = 1 is the user of interest. For convenience and 
simplicity it is assumed in the simulations that the propagation delays of the different users 
10 are multiples of Tc. Note that this assumption is not necessary if the sampling rate is higher 
than 1/Tc. 

First, the performance of the asynchronous CDMA receiver is found imder the 
assumption that the covariance matrix is known at the receiver. Even though an exact 
estimation of the covariance matrix is not possible to achieve practically in a real 

15 commimication system, this experiment provides an upper bound on the performance of the 
asynchronous CDMA receiver. The algorithm in Table 1 is used for this case. 

In FIG. 11, the probability of correct synchronization and the probability of correct 
demodulation are presented as a function of the signal-to-white noise ratio (Eb/No) under the 
perfect power-control scenario, i.e., all users have the same signal power (NFR = 0 dB). In 

20 this case the dimension of the full-rank receiver is equal to N = 31. Demodulation is 
accomplished by the use of the estimate of the propagation delay that is obtained firom the 
synchronization system. It is shown in FIG. 1 1 (NFR = 0 dB) that the detector performs 
better than the conventional detector. Actually, the conventional detector can be considered 
to be the rank 1 version of the proposed receiver. It is demonstrated in this figure that a rank 

25 2 filter already gains a significant improvement over the conventional detector. This 
improvement can be even more when the near-far ratio (NFR) is higher. A rank 3 filter 
achieves almost the same performance as the full rank (rank 31) filter. 

FIG. 12 demonstrates the probability of correct synchronization and the probability of 
correct demodulation in terms of Eb/No under the near-far situation of NFR = 3 dB. All of 

30 the other parameters are identical to those used in FIG. 11. Again the covariance matrix is 
assumed to be known. When the MAI becomes more severe, the receiver performs 



NFR = 




(55) 



A' 



2 ' 
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significantly better than the conventional detector. However, a greater rank (rank 4) is 
necessary in order to achieve an almost full-rank performance. 

The performance of the adaptive implementation of the receiver for the batch algorithm 
(FIR), given in Table. 3, is shown in FIGs. 13 to 16 with different scenarios for each figure. 
5 The different scenarios used in these simulation are summarized in Table 5. 

Table 5 : Exemplary scenarios for testing the adaptive multistage decomposition. 





N 


K 


L 


NFR 


FIG. 13 


31 


6 
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OdB 


FIG. 14 


31 
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3N 


3dB 


FIG. 15 


31 


6 


6N 


OdB 


FIG. 16 


31 


6 


6N 


3dB 



In FIGs. 13(a)-13(c), (NFR = 0 dB and L = 1/3A0 it is shown that rank 2 of receiver 

10 performs much better than the conventional detector. The case of rank 2 provides almost the 
same performance as the full-rank case while rank 3 outperforms the full-rank case. This is 
due to a lack of training samples in the full-rank case needed to obtain a good estimate of the 
covariance matrix. As previously stated, the reduced-rank version of the proposed detector 
requires a smaller number of samples to achieve synchronization and demodulation. FIG 

15 13(a) shows the probability of correct synchronization vs. signal-to-white noise ratio for k=6, 
N=31, and NFR=0 dB. FIG. 13(b) depicts the probability of correct demodulation vs. signal- 
to-white noise ratio for k=6, N=3 1 , and NFR=0 dB 

However, for the more severe MAI (NFR = 3 dB) a greater rank is needed in order to 
achieve the same or better performance than the full-rank case. This is shown in FIGs. 14 (a) 

20 and 14(b). By examining FIGs. 13 to 16 one can note that the full-rank receiver needs a 
larger amount of sample support to perform close to the known covariance matrix case while 
the reduced-rank receiver requires a smaller number of samples. The required computations 
to realize the adaptive full-rank receiver is sizeable due to the necessary matrix inversions, 
especially when the processing gain iV^of the CDMA system is larger XhanN- 31. 

25 FIGs. 15a and 15B show the batch algorithm in Table 3 with number of training samples 

L -6N. FIG. 15a shows a graph for the probability of correct synchronization vs. signal to 
white noise ratio for K=6, N=3 1 , and NFR is 0 dB. FIG. 1 5b shows a graph for the probability 
of correct demodulation vs. signal to white noise ratio for K=6, N=3\, and NFR is 0 dB. 
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FIGs. 16a and 16B show the Batch algorithm in Table 3 with number of training samples L 
=6K FIG. 16a shows a graph for the probability of correct synchronization vs. signal to 
white noise ratio for ^=6, N=3l, and NFR is 3 dB. FIG. 16b shows a graph for the 
probability of correct demodulation vs. signal to white noise ratio for Ar=6, N^3l, and NFR is 
5 3dB. 

Finally, the performance of the iterative or lER algorithm of the receiver in Table 4 is 
shown in FIGs. 17 and 18 with a = 1/3A'^ and FIGs, 19 and 20 with a = 1/6A^. Also for 
these two cases the near-far ratios are set respectively to 0 dB and 3 dB. FIG. 19a shows the 
probability of correct synchronization vs. signal to white noise ratio for K=6, N~3l, and NFR 
10 =0 dB, and FIG. 19a shows the probability of correct demodulation vs. signal to white noise 
ratio for K=6, N=^3\, and NFR = 0 dB. FIG. 20a shows the probability of correct 
synchronization vs. signal to white noise ratio for K=6, N=3l, and NFR = 3 dB, and FIG. 20b 
shows the probability of correct demodulation vs. signal to white noise ratio for K=^6, N-31, 
andNFR = 3dB. 

15 For the equal power case (NFR = 0 dB), the results show that the reduced-rank iterative 

version of the receiver outperforms the full-rank receiver, especially when the number of 
training samples is small. However, when the MAI become more severe (NFR = 3 dB), the 
performance of the reduced-rank iterative recesiver is degraded. A larger number of training 
samples is needed to obtain a better performance. When compared to the batch algorithm in 

20 Table 3 the recursive algorithm is computationally and physically much simpler, but does not 
perform as well as the batch algorithm. This is consistent with the results reported in [6]. To 
obtain a better performance, the iterative algorithm requires an equivalent number of samples 
as the full-rank case but with a much simpler computation. A more sophisticated iterative 
adaptive-weight filter algorithm can be used to obtain a better performance. However, in all 

25 cases the reduced-rank adaptive filter bank outperforms the conventional matched-filter 
detector. 

In summary, in this invention an adaptive near-far resistant receiver for an asynchronous 
CDMA system is developed without the use of prior synchronization. It does not require a 
training period or data-free observations. Only knowledge of the spreading sequence of the 
30 desired user is necessary. These are the requirements needed by the conventional receiver of 
an asynchronous CDMA system. Due to the near-far resistant property of the receiver, this 
new CDMA receiver does not require the accurate power control that is needed in a 
conventional receiver system. 
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It is shown that the receiver obtains a better performance than that of the conventional 
receiver in all simulations. Though an optimum adaptive receiver that uses a full-rank 
covariance matrix can suppress multiple-access interference that comes from the co-existing 
users, its practical implementation generally is not feasible due to a lack of the sample support 
5 that is needed to estimate the covariance matrix and the substantial computations needed to 
invert this matrix. 

To remedy this situation a lower-complexity adaptive receiver is developed that is based 
on the concept of the reduced-rank multistage decomposition [16]. With this method one can 
design a lower-complexity receiver without a significant loss in performance as compared 

10 with the full-rank system. Also the number of samples, needed to initially estimate the filter 
coefficients, is considerably decreased. This fact translates into an enhancement of the 
bandwidth efficiency of the underlying communication system or the potential performance 
gain when the size of the sample support for the adaptation of the filter coefficients is limited. 
Finally, an iterative version of the filter that utilizes the multistage decomposition is 

15 developed. This makes possible a receiver that can be realized as a bank of iterative adaptive 
filters that are multiplexed in time sequence. By this means, the coefficients of the filter are 
updated every time new data are received. The implementation of the iterative version 
requires a substantially smaller memory size than the earlier FIR versions. Although, the 
iterative version does not perform as well as the batch version, it is of practical importance 

20 because of its much simpler realization. The convergence rate of the iterative filters used for 
the filter bank can be improved by the used of higher order iterative filters in this DS-CDMA 
system. However the bank of iterative adaptive filters developed in this invention evidently 
results in a very fast, efficient and effective near-far resistant asynchronous CDMA receiver. 
As demonstrated above, using only the first 2 stages of the proposed receiver (the rank 2 

25 receiver) gains significant improvement in all cases over the conventional detector (which is 
equivalent to the first stage of the proposed receiver) while also the computational complexity 
of the rank 2 filter is only minimally greater than the conventional detector. 

Simulation results show that the receiver of the present invention can significantly 
outperform the conventional DS-CDMA receiver that uses a standard matched filter for 

30 acquisition. When compared to the MMSE type receiver, the receiver can achieve a similar 
performance level without the requirement of synchronization or known propagation delays. 
Some applications for the present invention include digital wireless networks including home 
networking, and any other wireless application that uses code modulation. 
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It will be recognized by those skilled in the art that various modifications may be made 
to the illustrated and other embodiments of the invention described above, without departing 
from the broad inventive scope thereof It will be understood therefore that the invention is 
not limited to the particular embodiments or arrangements disclosed, but is rather intended to 
cover any changes, adaptations or modifications which are within the scope and spirit of the 
invention. , 
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